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ABSTRACT 

Aims. We present molecular line observations of the southwestern part of the IC 348 young cluster, and we use them together with NIR and 
mm continuum data to determine the distribution of dense gas, search for molecular outflows, and analyze the ongoing star formation activity 
in the region. 

Methods. Our molecular line data consists of C'^'OCl-O) and N2H^(l-0) maps obtained with the FCRAO telescope at a resolution of about 
50" and CO(2-I) data obtained with the IRAM 30m telescope at a resolution of 1 1". 

Results. The dense gas southwest of IC 348 is concentrated in two groups of dense cores, each of them with a few solar masses of material and 
indications of CO depletion at high density. One of the core groups is actively forming stars, while the other seems starless. There is evidence 
for at least three bipolar molecular outflows in the region, two of them powered by previously identified Class sources while the other is 
powered by a still not well characterized low-luminosity object. The ongoing star formation activity is producing a small stellar subgroup in 
the cluster. Using the observed core characteristics and the star formation rate in the cluster, we propose that that similar episodes of stellar 
birth may have produced the subclustering seen in the halo of IC 348. 
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1. Introduction 

IC 348 is one of the most studied young clusters. Its stel- 
lar population has been observed at different wavelengths 
from the IR to the X rays (E^ada & Lada 1995; Luhman et alj 
1998; Herbi2 1998: Naiita et al.lfeood iMuench et al. 112001 
Preibisch & Zinneckeri i2004l) . and its surrounding gas. 



part of the Perseus molecular cloud, has been mapped 
with different resolutions in mm-w a ve lines and contin- 
uum teachiller&Cernich arol Il986t 'Hatchell et al. ' '2005^, 
Enoch etal. 2006; Ridg e et al. I I2OO6: Sun et al. 2006; 
Kirk et al. 2006). The cluster lies at a distance of 320 pc (e.g., 
Herbig 1119981) and consists of more than 300 stars distributed 



with a core-halo str ucture over a region 20 arcminutes in di- 
ameter jMuench et al. H2003I) . Superposed to the smooth dis- 
tribution of stars, there is a population of stellar groups ("sub- 
clu sters") each containin g 5-20 stars within a radius of 0.1-0.2 
pc jLada & Lada II 19951) . Most stars in IC 348 have formed at 
a close-to-constant rate over the last 2-3 Myr, although some 
cluster members ma y be significantly older (Herbig 1998; 
iLuhman et al.ll998l) . 

Star formation in IC 348 seems to have ceased toward 
the center but it continues at some level near its south- 
ern border, where the cluster meets the molecular cloud. 



IStrom et al. I ( ll974 identified an IR object to the SW of IC 
348 (IC348-IRS1 hereafter), and further observations in the 
visible and IR by Boulard et al. ( 1995) have resolved it into 
a bipolar nebul a likely due to an embedded star with a disk. 
iBachiller et al. I ill987l) identified several dense cores in its 
vicinity, and lMcCaughrean et al. ]([l994) found additional sign- 
posts of recent star formation in the form of shock-excited 
H2 emission. The most prominent of these H2 features is the 
HH21 1 flow, which originates from source HH21 1-MM and is 



inCOandSiO(Gueth&Guilloteau 1999 


Chandler& Richer 1 


2001; 


Palau et al. 


2006; Hirano et al. 


2006). Another promi- 
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ne nt H? feature h as been recently associated with an outflow 
bv lEisloff"el etal.l ( E003i) . who found an additional mm source 
that seems to power it. T he NIR and optical observation s of 
lEisloff^el et all te003') and lWalawender et al. I (12005". "2004 re- 
veal multiple HH objects in the region and suggest that addi- 
tional young stellar objects (YSOs) lie embedded in the dense 
gas. 

In this paper we present the results from a survey of the 
southwest region of IC 348 (IC348-SW hereafter) aimed to 
study its star formation activity and the relation between the 
dense gas and the molecular outflows. These data reveal the 
presence of several aggregates of dense cores, some of them 
starless, together with at least three molecular outflows (two 
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newly mapped in CO), and help clarify the kinematics of both 
the outflow gas and the dense material. Combining our obser- 
vations with published IR data, we show that the subclusters in 
the halo of IC 348 first identified by Lada & Lada ( 1995) could 
have originated from star formation episodes inside small core 
aggregates like the one currently active in IC348-SW. 

2. Observations 

We observed the IC348-SW region in C"^O(l-0) and N2H+(1- 
0) with the (then) 16-pixel SEQUOIA array receiver on the 
FCRAO' telescope in 2001 April. The observations were done 
in frequency switching mode using a correlator that provided 
velocity resolutions of 0.03 km s"' (C'^O) and 0.06 km s"' 
(N2H^). The pointing was checked and corrected using obser- 
vations of SiO masers, and the data were converted into the 
main-beam temperature scale using an efficiency of 0.55. The 
FWHM of the telescope beam was approximately 50". 

We observed the vicinity of IC348-IRS 1 in the 1 .2mm con- 
tinuum with the MAMBOl bolometer array at the IRAM 30m 
telescope in 1999 December. One on-the-fly map was made 
with a scanning speed of 4" s a wobbler period of 0.1 s, 
and a throw of 41". The atmospheric optical depth was esti- 
mated from sky dips carried out immediately before and af- 
ter the map, and the data were reduced with the NIC soft- 
ware package. We also observed the IC348-SW region simul- 
taneously in CO(l-O) and CO(2-l) with the IRAM 30m tele- 
scope in 2000 October and 2001 October. The observations 
were done in dual-polarization, position switching mode and 
were centered at 0-2000=03:43:58.8, (J2000=32:01:51. The refer- 
ence position was 20 arcminutes north from the map center, 
and it was checked to be free from detectable emission in the 
velocity range of interest. A correlator split into four sections 
provided velocity resolutions of 0.20 and 0.41 km s"^ for the 1- 
and 2-1 spectra, respectively. The atmosphere was calibrated 
by observing ambient and cold loads, and standard efficiency 
values were used to convert the telescope intensities into the 
main-beam brightness scale. The pointing was checked mak- 
ing cross scans on bright continuum sources, and the FWHM 
of the telescope beam was 21" at the CO(l-O) frequency and 
11" at the CO(2-l) and 1.2mm continuum frequencies. 

Near-infrared observations were made using the 3.8m 
UKIRT with the UIST aiTay camera on the night of 2002 
December 5 under excellent seeing conditions (0.5" in K 
band). The camera provided a plate scale of 0.12" per pixel 
with a total field of view of 120" x 120." We obtained a 9 point 
jitter of 60 second exposure time in the K band. Standard proce- 
dures for data acquisition and reduction were followed, involv- 
ing flat fielding, sky, and dark subtraction of the raw frames. 
The limiting magnitude of the image was 17.6. 

3. Large-scale distribution and dense cores 

Figure 1 shows in contours our large-scale FCRAO maps 
of the C'**O(l-0) (left) and N2H+(l-0) (right) emission to- 

' FCRAO is supported in part by the National Science Foundation 
under grant AST 94-20159, and is operated with permission of the 
Metropolitan District Commission, Commonwealth of Massachusetts. 



ward the south-west vicinity of the IC 348 cluster super- 
posed on th e 2MASS K-band and Spitzer IRACl (3.6 //m) 
images (see lj0rgensen et al. I l2006l for full IRAC maps of 
Perseus). As the C^^O map shows, most of the molecular 
material lies along an arc whose center approximately coin- 
cides with the center of the IC 348 cluster and has a radius 
of about 12' or 1 pc (see flatchefl et al. 200^ for a com- 
plete C'^O(l-O) map of the region). The dense gas, traced 
by N2H^(l-0) in the right panel, consists of two groups of 
cores along the C'^0 arc plus a weaker core to the north, 
close to the IC 348 center The western group of cores, which 
we will refer to as IC348-SW1, was mapped previously in 
NH3 by Bachiller et al. ( 1987). It has an approximate horse- 
shoe shape and consists of three cores named A, B, and C by 
iBachiller et al. 1 lll987l) (see Fig. I). Some of these cores con- 
tain well known young st ellar objects ^(YSO s): core B is asso- 
ciated with IC 348-IRS1 ( Strom etal. '1974 *), core C contains 
HH21 1-MM (nMcCau ghrean e t al.iiil994: iGueth & Guilloteaul 
fl999), an d between cores B and C lies the mm source IC348 
MMS of jEisloff'el et al. II2003I) . The second group of cores is 
located to the southeast of IC348-SW1 and will be referred to 
as IC348-SW2. It consists of two cores, none of them associ- 
ated with a known YSO or an IRAS source. As Fig. 1 shows, 
both IC348-SW1 and SW2, and their suiTounding molecular 
gas, coincide with an almost total absence of scattered light at 
3.6 jum. This suggests that the two regions lie in the front part 
of the IC 348 cluster. 

The N2H^ emission of IC348-SW1 in Fig. 1 is remarkably 
similar to the NH3 emission mapped by Bachiller et al . llll987B 
and to the continuum emis sion rti a pped by .Hatchell et al. I 
(|200j at 850 jum and Enoc h et al. I ilOOdti at 1.1 mm. Such 
a good agreement between maps suggests that these tracers 
reflect the true distribution of dense gas in the region. The 
C'^O(l-O) emission, on the other hand, shows little contrast 
over the map and and peaks at a different position. The bright- 
est C'**O(l-0) is located between the A and B cores of IC348- 
SWl, and the IC348-SW2 group of N2H+ cores coincides wifli 
a region of relatively weak C'**0 emission. This contrast be- 
tween the C'^O maps and the maps of N2H+, NH3, and the 
continuum is unlikely to result from optical depth effects, as 
the C'**O(l-0) lines are less than 4 K in intensity and there- 
fore not optically thick (the thick CO(2-l) lines are typically 
20 K bright). It most likely results from C"*0 being strongly 
depleted at the highest densities. C'^O (and CO) depletion due 
to freeze out onto dust grains is a common phenomenon in the 
cold, low-turbulence cores o f clouds Uke Taurus , were it is by 
now well characterized C^g.- ICaselli et al. Jl999; Tafa lla et al. I 
l2002h . Its finding in the warmer and more turbulent cores of 
IC348-SW indicates that CO depletion occurs at a relatively 
broad range of conditions. 

To estimate the mass of each dense core we use the N2H^ 
emission. We first estimate the central H2 column density 
NiHi) by assuming an N2H^ excitation similar to the one found 
in Taurus, and an N2H^ abundance of 1 .5 x 10 (also as found 
in Taurus, e.g., Tafalla et al. 2004). This method produces H2 
column densities for cores A, B, and C that agree wit hin a fac- 
tor of 2 with the continuum-based estimates by Enoch et al. I 
(■2006t) . and we take this factor as indication of the level of un- 
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Fig. 1. FCRAO maps of the molecular gas in the south-western part of the IC 348 cluster. Left: C"*O(l-0) integrated intensity 
map in contours superposed on the K-band 2MASS image. The C"*0 emission forms an arc that partly surrounds the IC 348 
cluster (located at the top left). Right: N2H^(l-0) emission (contours) superposed on the Spitzer IRACl (3.6 /im) image. The 
dense gas traced by NiH^ forms two groups of dense cores, IC348-SW1 and IC348-SW2, and a small condensation toward 
the cluster center The map center is at Qf200{)=03:43:58.8, (52000=32:01:51, and the first contour and interval is 0.5 K km s~^ 
for C'^O(l-O) and 1 K km s"' for N2H^(l-0). The dashed line indicates the region covered by the FCRAO observations. The 
C'**O(l-0) map has been convolved to a final resolution of 65" (FWHM) to filter out high frequency noise. 



certainty of our estimate. We then measure the core radius R 
from the N2H+ map and assume an i nternal density stru cture of 
a critical Bonnor-Ebert sphere (e.g., lAlves et al. I2OO1I) . In this 
way, the mass of the core is given by 

(i^Kofi^f^- 

With this approximation, we derive masses of 9, 7, and 3 
Mq for cores A, B, and C, respectively, which are in rea- 
sonable agreement w ith the NHs-based mass estimates of 
iBachiller etall ( Il987l) . For the IC348-SW2 group of cores, 
we derive masses of 5 (western core) and 1 .5 (eastern 
core), and for the core close to the center of the IC 348 cluster, 
we estimate an approximate mass of 1 M©. 

We have explored the velocity structure of IC348-SW1 and 
SW2 using the N2H^(l-0) data, and we illustrate the main 
features with the two position-velocity (PV) diagrams of Fig. 
2. The PV diagram for SWl follows a curved path along the 
horseshoe, and the diagram for SW2 follows the long axis of 
the region (PA = 60°). As the top panel shows, the A core 
presents a single velocity component with a small velocity gra- 
dient («i 1.3 km s ' pc ') and a peak at around 8.6 km s The 



B and C cores, on the other hand, present lines with two peaks 
separated by about 0.5 km s"', being brighter the blue compo- 
nent in core B and the red component in core C. The presence 
of these two components seems not to result from an overlap 
between the cores, as each core appears distinct and centrally 
concentrated in the map of Fig. 1 . In addition, the weaker red 
component of core B peaks at the same position as the brighter 
blue component, as if the two were correlated, and a similar 
but weaker correlation can be seen in core C. This behavior 
suggests that in each core the two components have a common 
origin. As we will see in the next section, the B and C cores are 
affected by outflows that have akeady accelerated part of the 
ambient gas (as seen in CO). It is therefore likely that the N2H^ 
components arise from outflow acceleration, as it has been pre- 
viously found in other systems like L1228 (iTafalla & Mvers I 
■ 1997i) . 

The velocity structure of SW2, on the other hand, is simpler 
than that of SWl. The N2H^(l-0) line presents a single peak 
over the region, and there is a velocity gradient again at the 
level of 1 .3 km s"' pc"' . This velocity gradient could represent 
a smooth change in the central velocity over SW2 or result from 
the two cores of SW2 having different velocities at the level of 
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Fig. 2. Position-velocity diagrams of the N2H+(JFiF= 101-012) 
("isolated component") emission toward IC348-SW1 (top) and 
IC348-SW2 (bottom). For IC348-SW1, the position axis fol- 
lows a curved path that intersects the peaks of the A, B, and 
C cores, and the the approximate location of each core is in- 
dicated with arrows to the right. For IC348-SW2, the position 
axis follows a straight line with PA - 60°. In both panels, the 
first contour and interval are 0. 1 K. 

0.2 km s"'. The contrast with the double-component velocity 
structure of SWl may result from the SW2 region having no 
evidence for star formation and or outflows (see below). 

4. Molecular outflows 

We have searched for molecular outflows the IC348-SW1 and 
IC348-SW2 regions using CO(l-O) and CO(2-l) observations. 
No outflows were found towards IC348-SW2, while at least 
three outflows (two newly detected in CO) were found in 
IC348-SW1. To simplify the presentation, we separate the out- 
flow velocity regime into "fast" and "slow." 

4.1. Fast gas 

We define as "fast" the gas that moves at more than 5 km s"' 
with respect to the ambient cloud (ambient Vlsr = 8.6 km s"' 
in IC348-SW1). This 5 km s"' value corresponds to ten times 
flie typical FWHM of the N2H+ line in IC348-SW1, and there- 
fore guarantees the absence of contamination by ambient emis- 
sion. The exact choice of the value, however, has little efi'ect on 
the following discussion. 

As Fig. 3 shows, the fast regime is dominated by two 
bipolar outflows, each centered on a bright 1.2mm compact 
source. The southern o utflow is the well-known HH 211 sys- 
tem, first identified by iMcCaughrean et a~ ] (I1994 from its 
H2 emission , and later mapped at high resolution in CO 
and SiO bv iGueth & Guilloteau I (Il999t. IChandler & Richerl 
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Fig. 3. Fast outflows in IC348-SW1. The contours show the 
CO(2-l) emission between 5 and 25 km s"' with respect to 
ambient cloud (at Vlsr = 8.6 km s"'). The solid blue contours 
represent blueshifted gas and the dashed red contours represent 
redshifted gas. The grey scale shows the 1.2 mm continuum 
emission. The two bright mm sources at the intersections of 
the red and blue contours are IC348-SMM2 (top) and HH21 1- 
MM (bottom). Offsets as in Fig. 1, and first contour and contour 
interval is 3 K km s"' . 



pea k at its center have been v ariously refer red as HH211-mm 
by d McCaughrean et al. lfl994li) or SMMl (iWalawender et al. I 
2006). The second fast outflow coincides w i th the extended 
H2 emission detected by^ McCaughrean et al. 1 (ll994h (flieirH2 
"chain") and lEisloff^el et al. ■ (,2003J (their "region 1"), and it 
also coincides with t he HH797 A and B featu res recently 
observed in [Sii] by IWalawender et al. ] i2005). Compared 
with the optical and IR data, our CO observations pro- 
vide velocity information, and show that the outflow is 
bipolar with respec t to the north ern mm source (referred 
as IC348 MMS bv lEisloffel et al. . .2003. and as SMM2 by 
IWalawender et al. II2006I) . The CO data, therefore, shows un- 
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of the HH211 outflow is unresolved with arcsecond beams. 
The SMM2 outflow, on the other hand, is unresolved by our 
11" single-dish beam at the highest speeds. In addition, its 
CO spectra present at some positions secondary peaks of 
high velocity (see Fig. 4). These peaks, often called "bul- 
lets," are commonly associated with jet-like components in out- 
flows (Bachiller et al. 1990), and their presence in the SMM2 
system suggest the existence of extremely collimated gas. 
Interferometric observations of the SMM2 outflow are neces- 
sary to further characterize this component. 
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Fig. 4. CO(2-l) spectra and position velocity (PV) diagram 
along the IC348-SMM2 outflow. The spectra present secondary 
peaks of high velocity emission ("bullets") at positions sym- 
metrically located with respect to SMM2. The PV diagram 
shows how the outflow speed increases with distance from 
SMM2 and how the bullets are highly localized in both veloc- 
ity and position. Note that the offset in the spectra are referred 
to the map center (Fig. 1), so SMM2 is at (-25", 75"), while 
in the PV diagram, the position axis measures distance from 
SMM2. The dashed lines mark the boundary of fast regime. 
First contour an interval in the PV diagram are 0. 1 K. 

am biguously that SMM2 po wers the outflow (we will follow 
flie IWalawender et al.1l2006l notation in the following discus- 
sion). 

The HH211 and SMM2 outflows present several differ- 
ences and similarities. HH211 is more compact than the 
SMM2 outflow, and although this may indicate an intrinsic 
difference between the outflows it could also result from a 
projection eff'ect. Indeed, the SMM2 outflow presents some 
mixing of blue and red emission at low velocities (see be- 
low), which suggests that the flow runs close to the plane 
of the sky. In any case, both outflows seem to belong to 
th e class of highly coll imated flows. The interferometric data 
oflGueth & Gu iUoteau (.1999.) (also.Chandler & Richer 1200 it 
IPalau et al. 12006: Hirano et al. l2006l) show that the fastest part 



4.2. Slow gas 

The threshold of 5 km s ' used for the fast gas guarantees the 
absence of contamination by low-velocity material, but limits 
our sensitivity to weak or slow outflows. As a second step in 
our outflow search, we now explore the gas moving at veloci- 
ties lower than this threshold, but still avoiding contamination 
from ambient gas. We have experimented with different veloc- 
ity ranges and have found that the gas moving between 3 and 
5 km s"' from the ambient cloud still seems free from ambient 
contamination. For this reason, we will refer to this range as 
the "slow gas" regime, although it should be noted that this gas 
still moves supersonically with respect to the ambient cloud. 

Figure 5 shows the slow gas regime of CO(2-l) toward the 
IC348-SW1 region (solid blue and dashed red contours) super- 
posed to the N2H^(l-0) emission (dotted lines) and the Spitzer 
Space Telescope 24 ^m image from the Cores to Disks (c2d) 
Legacy Project dEvans et al. Il2003t) . The HH211 and IC348- 
SMM2 outflows are still the dominant elements of the slow 
CO emission, but a number of new features are visible. In the 
IC348-SMM2 outflow, there is a red component at the end of 
the northern blue lobe and a blue arc to its northeast. The ori- 
gin of these components is suggested by their location with 
respect to the dense gas traced by N2H^. The N2H^ emis- 
sion shows that the SMM2 outflow runs between the A and 
B cores and that the two anomalous features occur at the edge 
of the region, where the outflow leaves the dense gas and en- 
counters the surrounding medium. The anomalous red emis- 
sion continues the direction of the flow and coincides with 
a region where the sho cked H? emiss ion becomes b righter 
and broade r (McCaughr ean et al 1119941 ^loff'el et aLl boOSt 
IWalawen'der et al. 12 006). This suggests that the anomalous red 
CO emission results from a broadening of the outflow when 
it encounters the surrounding lower density medium, as it has 
been o bserved in other systems like the Orion HH212 outflow 
jLee et al. ■■2000.) . The outflow broadening may result from a 
change in the gas regime from isothermal to adiabatic, and if 
the outflow moves close to the plane of the sky, it can naturally 
produce the obser ved superposition of blue and red gas (e.g., 
ICabrit et al. Il988l) . The blue arc to the east, on the other hand, 
has a less clear origin. It could result from a strong deflection 
of the outflow, but there is no evidence for an obstacle along 
the outflow path. In fact, the wind responsible for the outflow 
seems to continue unimpeded past the region with anomalous 
blue emission. This is inferred from the long and highly colli- 
mated chain of H2 emission, which continues at least up to the 
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Fig. 5. Slow outflow gas toward IC348-SW1. The contours show the CO(2-l) emission integrated between 3 and 5 km s"' with 
respect to the ambient cloud (at Vlsr = 8.6 km s '). The solid blue lines represent blueshifted emission, and the dashed red 
lines represent redshifted emission. The black dotted contours indicate the N2H^(l-0) emission also shown in Fig. 1, and the 
grey-scale image in the background is the 24 yum Spitzer MIPS image from the c2d Legacy Project. At least three CO bipolar 
outflows are present in the region (see text). Off'sets as in Fig. 1, and first CO contour and contour interval is 1 K km s"^. 



bright "1-w" knot of lEisloffel et al. I ( l2003h . Such a knot is lo- 
cated at offset (-54", 243") with respect to our map center, and 
therefore is well aligned with the collimated CO outflow past 
the region of anomalous blue gas (it is in fact so much further 
north from the blue gas that it lies outside our map). A more 
likely interpretation, therefore, is that the blue gas is not re- 
lated to the outflow. Fig. 5 shows that the blue gas forms an arc 
parallel to the dense gas distribution traced by NaH^ and that 
approximately points towards the center of the IC 348 cluster 
(about 1 pc away in projection). We have seen from the C'^O 
emission (Fig. 1) that the cluster seems to have excavated a cir- 
cular hole in the molecular gas, so it is possible that the blue 
CO emission results from this destructive action of the cluster 
and not from the CO flow. This would also explain naturally the 
shift to the blue, as the IC348-SW1 region is located slightly in 
front of the cluster (Section 3), so the dense gas will be pushed 
from behind. 

After the CO emission associated directly with the HH21 1 
and SMM2 outflows, the next brightest peak of the slow regime 
occurs near (30", -40") and is part of a region of red gas 
that lies along the southwestern edge of core B (Fig. 5). 



This red gas is not continuously connected to the red lobe of 
SMM2, but the presence of another region of red emission near 
( -20", -20") and a number of nearby H2 knots in the images 
of lWalawender et al. Ill200 6^ suggests that the two pieces of red 
gas are indeed connected and associated with the southern lobe 
of the SMM2 outflow. If so, the southern lobe would be ap- 
proximately as long as the northern lobe. 

About 40" north of the bright red peak just discussed, 
there is a triangular patch of blue gas whose northern ver- 
tex coincides with the "region 3" of H2 emission in the study 
of .Eisloffel et al. (20(13). As Fig. 5 shows, this region also 
presents bright 24 jum emission, which suggests the presence 
of an additional embedded object. 

A better defined feature of the slow CO emission occurs 
near (-100", 100"). There, a blue and a red lobe of gas form 
a very compact outflow at the eastern edge of core A. This 
newly found CO flow runs approximately north-south, paral- 
lel to the SMM2 outflow, and seems associated with t he "re - 
gion 2" of H2 emission identified by teisloff'el et al. I ll2003h . 
The close-up view in Fig. 6 shows that the outflow has dif- 
fuse H2 emission both towards its north and south lobes, and 
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Table 1. Outflow energetics 




Fig. 6. IC348-SMM3 outflow. The solid (blue) and dashed (red) 
contours indicate the CO(2-l) emission blueshifted and red- 
shifted between 3 and 5 km s ' with respect to the ambient 
cloud (slow outflow regime). The grey scale shows our UKIRT 
K-band image and the plus sign indicates the position of a 24 
yum source. Offsets as in Fig. 1, and first contour and interval is 
l.OKkms-i. 



Outflow 


Mass 


Momentum 


Energy 




(Mo) 


(Mokms-') 


(erg) 


HH211-MM 


1.5 10-2 


1.5 10-' 


6 lO''^ 


IC348-SMM2 


3.3 10-2 


3.1 10' 


1 10** 


IC348-SMM3 


10-3 


4 10-3 


4 10*1 



7), and these two lines fall inside the 4.7 yum-wide bandpass of 
the 24 jum MIPS detecto r. Bright [Si] emission at 25.2 /im has 
in fact been observed bv lNoriega-Crespo et al. I (l2004l) toward 
the Cepheus E outflow, which is also powered by a very young 
object, so this emission is also expected from the outflows in 
the IC348-SW1 region. Until spectroscopic observations of the 
24 yum sources are available, their association with a YSO re- 
mains uncertain (the lower wavelength Spitze r bands of IRAC 
can also be contaminated by H2 emission, |j0rgensen et al. I 
[2006). A more secure counterpart for the exciting source of the 
compact outflow is provided by the submm c ontinuum emis- 
sion observed by different authors. Source 15 of lHatchell et al. I 
(2005), source 102 of Enoch et al. (20(^3), and source SMM3 
of lWalawender et al. I ( l200dl) all coincide with the approximate 
center of the CO outflow, and even our 1 .2 mm observations, 
where the source is close to the map edge, show evidence of 
emission at the 100 mJy level. Given this more secure identi- 
fication of the source at submm wavelengths, from now on we 
will refer to the outflow driving source as IC348-SMM3. 



that its center lies about 8" to the northeast of 2MASS source 
03435056+3203180. The 2MASS source, however, is probably 
not the outflow source, as its JHK co lors are cons istent with 
those of a background or TTauri star dLada & Adams ■■1992) . 
A more likely candidate for the exciting source is provided by 
the lAjim Spitzer image. As Fig. 5 shows, there are two sources 
superposed to the CO emission, one toward the geometric cen- 
ter of the outflow (plus sign in Fig. 6) and the other coinciding 
with the bright H2 knot at the end of the red lobe. Both 24 yum 
sources are different from the 2MASS source and are classi- 
fied a s YSO candidates in the current c2d catalog (Evans et al. 
I2OO5I) . The northern source, close to the outflow center, is most 
likely the driving source of the outflow, while the second source 
may result not from an embedded object but from line emis- 
sion in the shock. A similar situation seems to occur toward the 
HH2 1 1 outflow, where the end of the blue lobe coincides with a 
bright 24 jim peak (Fig. 5). This 24 fim peak is also classified as 
a YSO candidate in the c2d catalogue, but it has no evidence for 
mm or submm emission (e.g.. Fig. 3), while its position coin- 
cides with the p oint of strongest IR and o ptical shock emission 
in HH211 (^g.- IWalawender et al. IE0O6I) . Line emission from 
outflow shocks near 24 iim is in fact expected from theoreti- 
cal grounds.flollenbach & McKee ( 1989) predict that the [S i] 
line at 25.2 jim and the [Feii] line at 26.0 /im will be among 
the brightest atomic fine-structure lines in shocks at densities 
between 10"^ and 10^ cm-^, typical of dense cores (their Fig. 



4.3. Outflow energetics and nature of the powering 
sources 

We calculate the mass, momentum, and energy of the outflows 
using the CO(2-l) emission, which we assume is optically thin 
in the outflow regime and that has an excitation temperature 
of 20 K, as suggested by the ambient peak intensities (these 
two assumptions make our estimate a lower limit). We als o as- 
sume a CO abundance of 8.5 10-^ Jprerking et al. lfl982i) and 
make no correction for outflow inclination. Neglecting the con- 
tribution of gas moving at less than 3 km s"' from the ambient 
cloud (Vlsr - 8.6 km s '), we obtain the outflow parameters 
listed in Table 1 . Although the HH2 1 1 outflow has been stud- 
ied before, previous estimates of its parameters were based on 
interferometer data, and therefore missed a significant fraction 
of the emission. The Table 1 values, therefore, contain the first 
full estimate of the HH21 1 energetics. 

As can be seen in the table, the HH211 and SMM2 out- 
flows have very similar parameters. Their mass, momentum, 
and energy differ by a factor of 2 or less, which is little given 
all the uncertainties involved. The slightly larger values of 
the SMM2 flow (despite its location closer to the plane of 
the sky) may arise from the larger luminosity of the source: 
8 Ln versus the 4 L n of HH211-MM (Froebric h'etaLll2003[ 
iRisloffel et al"1 l2003V These similar energetics, together with 
the similar speeds and high collimation, suggest that the pow- 
ering sources of the HH211 and SMM2 outflows share simi- 
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lar phy sical p roperties and evolutionary state. iFroebrich et al. I 
(l2003h and Eislof fel et al. I (l20 03l) have classified these sources 
as Class (Andre et al. Ill99 3). and the outflow parameters in 
Table 1 confirm that the two objects belong to the youngest 
protostellar phase. 

The IC348-SMM3 outflow, on the other hand, has a very 
small energy content, close to the detection limit of our survey. 
Compared with the HH211 and SMM2 outflows, the SMM3 
outflow is about one order of magnitude less massive and has 
two orders of magnitude less kinetic energy. These numbers, 
together with the small spatial extent of the CO emission (about 
30" or 0.04 pc) suggest that the flow is powered by a source of 
much lower luminosity than HH21 1-MM or SMM2. The exact 
luminosity of SMM3, unfortunately, is not well constrained. 
The 24 fim flux may be contaminated by outflow emission (see 
above), and no data close to the peak of the SED are available 
yet (the IRAS data are confused by the extended emission from 
the IC 348 cluster). Although the outflow parameters suggest 
that SMM3 is a low luminosity source, it is unlikely that its 
energy output is as low as that of the very low lu minosity object 
(VELLO) L1014-IRS (0.09 L^. lYoung et'an[2004) . For this 
object, ^ourke et al. (2005) estimate an outflow mass about 
70 times smaller and a kinetic energy at least two orders of 
magnitude lower than those of the SMM3 outflow. Although 
the interferometric data used to derive the L1014-IRS outflow 
parameters only provide a lower limit because of filtering of 
the extended emission, the difference in values is too large to 
be explained as an artifact. It seems therefore likely that SMM3 
is more luminous than 0.09 Lq. 

5. Star formation in IC348-SW and its relation to 
the IC 348 cluster 

We have seen that the IC348-SW1 group of cores is relatively 
rich in young steflar objects. Sources HH211-MM, IC348- 
SMM2, and IC348-SMM3 all drive bipolar outflows and are 
therefore YSOs likely formed in IC348-SW1 wit hin the last 
10^ yr (Class objects are even younger, lAndre et al. I 
The IC348-IRS 1 object, whose discovery provided the 
first indication that star formation was still ongoing in IC 348 
JStrom et al. Ill974l) . has a less clear status. It has been inter- 
preted as a late B star with an edge-on disk b vlBoulard et al. I 
(I1995I) . but recent IR observations by Luhm an et al. 



suggest that it has an early M spectral type. Its l arge extinction 
(iLuh man et al. I I 199^ and reflection nebulositv lBoulard et al. I 
(Il995) indicate youth, but the lack of a well defined outflow 
(Section 4.2) suggests that it is more evolved than the other 
three sources; a complete spectral energy distribution is still 
needed to properly classify this YSO. Other y oung objects in 
the region have bee n ident ified in the surveys of Herbig ( 1998) 
and iLuhman et al" 1 (ll998l) . who searched for evidence of H„ 
emission and IR excess. Among these objects, IfAHA 7 and 
8 stand out for lying projected onto the N2H^ emission and 
having Hq, equivalent widths of classica l T Tauri stars (35 and 
30 A for IfAHA 7 and 8, r espectively. iHerbig Ill998l) . Given 
their expected age («; 10^ vr. lLuhman et al. Il998l) . however, it 
is unclear whether these stars have originated from the IC348- 
SWl group of cores or they belong to a more distributed stellar 



population that has already drifted from its natal place. Even 
discounting these older objects, the presence of 4 YSOs, two 
of them of Class 0, makes IC348-SW1 the most active star- 
forming site of the IC 348 cluster 

The study of IC348-SW1 may h elp elucidate the proc ess of 
star formation in the IC 348 cluster. iLada & Lada I d 19951) have 
found that the stellar distribution in the outer cluster presents 
significant substructure. In addition to a central core of stars 
with a radius of 0.5 pc, the authors identify eight subclusters, 
each containing 5-20 sources within a radius of 0.1-0.2 pc. 
The sizes of these subclusters match the approximate size of 
the IC348-SW1 and SW2 regions as measured from the N2H'^ 
maps (radii 0. 1-0.3 pc), and this suggests a connection between 
the two types of objects. Using the 4 YSOs identified in IC348- 
SWl, which lie inside a circle of 85" radius, we crudely es- 
timate a stellar density of 70 stars pc"^. This density is al- 
ready half the mean stellar density of the Lada & Lada subclus- 
ters, and equals the density of their subcluster "c." Adding the 
nearby sources IfAH A 7 and 8 (Herbig 1998) and sources 49, 
160, and 124 from the lLuhm'anI d 1999ft catalog further enriches 
the stellar census toward IC348-SW1, to which new sources 
will be incorporated if the cores form additional stars with the 
dense gas still available. The final result of star formation in 
IC348-SW1 and its vicinity will therefore be a stellar density 
enhancement similar to the subclusters found bv lLada & Ladal 
(I1995I) . 

To explore whether star formation in regions like IC348- 
SWl is a viable mechanism to produce the observed subcluster- 
ing, we need to investigate whether the number of stars found 
in subclusters is consistent with their production rate and dis- 
persal time scale. A given group of stars will appear as a dis- 
tinct density enhancement until its members drift apart due to 
proper motions, and for the case of IC348-SW1 we can es- 
timate this timescale from the typical velocity difference be- 
tween its constituent cores (the linewidths may be contami- 
nated by core disruption, see section 3). From our N2H^(1- 
0) spectra, we estimate that the typical core-to-core velocity 
diff'erence i s less than 0.3 km s~' in agreement with the NH3 
estimate of IBachiller et al. I (Il987l) . This line-of-sight estimate 
implies a plane of the sky dispersion close to 0.4 km s ' if 
the velocity field is isotropic. As the su bclusters have a typical 
radius of 0.15 pc iLada & Lada I ( 1 19951) . stars will travel a di- 
ameter distance in about 0.75 Myr, which we take as a typical 
stellar group lifetime. To calculate the number of stars formed 
in this time, we assume that the cluster has been formi ng star s 
at a constant rate over the last 3 Myr (Muench et ain i2003l) . 
and that half of its 380 stars have formed in the outer clus- 
ter (Lada & Ladalll995l) . Thus, over the last 0.75 Myr, about 
47 stars have formed in the outer IC 348 cluste r. This number 
should be compared with the 82 stars found bv iLada & Lada I 
( 1995) in the subclusters of the halo, but first we need to sub- 
tract the contamination by the more diffuse cluster population. 
Dividing the number of non subcluster objects in the outer clus- 
ter by the area of this region, we estimate that the contamina- 
tion fraction is 1/3 (density of diffuse population is 50 stars pc"^ 
and mean density in a subcluster is 150 stars pc"^). This implies 
that there are about 55 stars in the subclusters of the outer IC 
348, which is roughly consistent with the 47 stars expected if 
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the stars formed in regions like IC348-SW1. The subcluster- 
ing structure, therefore, could have resulted from star forma- 
tion in aggregates of dense cores like the ones we have studied. 
Indeed, if star formation in these loose aggregates produces a 
different proportion of low mass stars than the formation in a 
more compact environment, likely responsible for the central 
stellar core of IC 348, the difference may help expl ain the ra- 
dial dependence of the IMF found bv iMuench et al. il(2003l) . 

Finally, we note that the velocity dispersion between aggre- 
gates of cores seems to be larger than the velocity dispersion 
between the cores of a given aggregate. Our N2H^ mapping of 
the IC 348 vicinity is not complete enough to determine the 
full statistics of the aggregate-to-aggregate kinematics, but it 
already provides some interesting hints. The SWl, SW2 ag- 
gregates and the core towards the IC 348 center have mean 
LSR velocities of 8.6, 9.1, and 7.9 km s ^ respectively, from 
which we estimate a (necessarily crude) rms value of 0.6 km 
s"' together with a trend for the velocity difference to increase 
with distance between aggregates. This dispersion is larger than 
the difference of about 0.3 km s~' found between the cores of 
IC348-SW1, and suggests that most of the initial velocity dis- 
persion between the stars in the core will result from the rel- 
ative velocities of the aggregates and not from the small scale 
turbulence inside the aggregates. It also suggests that the stars 
formed from a given aggregate may keep a distinct velocity pat- 
tern for about 0.75 Myr, which should be testable with a proper 
motion study of the subclusters. 

6. Summary 

We have mapped the molecular gas south-west of the IC 
348 young cluster in Ci^O(l-O), N2H+(l-0), and CO(2-l). 
Combining these data with observations of the young stars and 
YSOs in the NIR and mm continuum, we have reached the fol- 
lowing conclusions. 

1 . The dense gas south-west of IC 348 is concentrated in 
two regions, referred here as IC348-SW1 and IC348-SW2, 
each containing several cores of a few solar masses of mate- 
rial. The two regions are embedded is lower density gas and 
they seem located in the front part of the IC 348 cluster. While 
IC348-SW1 is actively forming stars, IC348-SW2 seems qui- 
escent. 

2. The distribution of N2H^(l-0) emission is similar to the 
distribution of NH3 and submm continuum mapped by previ- 
ous authors, but differs significantly from that of C"^O(l-0). As 
the rare isotope is not optically thick, the disagreement suggests 
that the CO molecule is depleted at the core centers (likely due 
to freeze out onto grains) as previously seen in the more quies- 
cent cores of Taurus. 

3. The fast CO emission (|Ay| > 5 km s"' with respect to 
the ambient cloud) is dominated by two highly coUimated out- 
flows, each associated with a bright mm/submm source. The 
first outflow is the molecular counterpart of HH211 and has 
already been mapped using interferometers by a number of au- 
thors. The second outflow, not previously mapped in CO, is 
associated with a previously known chain of H2 emission and 
several HH objects, and is referred here as the IC348-SMM2 
outflow. The characteristics and energetic content of the two 



outflows are consistent with them being excited by Class 
sources. 

4. At low velocities (3 < |Ay| < 5 km s"'), the CO emission 
presents additional features associated with the SMM2 outflow 
and (probably) the interaction of the cluster and the dense gas. 
In addition, we find a compact outflow associated with with a 
mm/submm object (SMM3) also detected at 24 //m. The small 
size and low energy content of this new outflow suggests that it 
is powered by a source of rather low luminosity. 

5. The ongoing star formation activity in IC348-SW1 is 
producing a significant stellar density enhancement, similar to 
others previously identified in the halo of the IC 348 cluster 
from NIR observations. From an estimate of the production rate 
and dispersal timescale of aggregates like IC348-SW1, we con- 
clude that similar episodes of star formation could have origi- 
nated the substructure seen in the IC 348 cluster 
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